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ABSTRACT 


A  new  method  is  proposed  for  improving  thfj  body-wave 
magnitude  determination  by  using  the  observed  values  of 
the  body-wave  magnitude  (m^)  together  with  the  first  motion 
directions,  to  ootain  by  least  squares  analysis  the  best 
double  couple  source  parameters;  the  resulting  radiation 
pattern  is  then  integrated  spatially  to  provide  a  corrected 
estimate  of  the  magnitude.  Results  for  a  number  of  events 
previously  studied  by  other  investigators  are  presented. 
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INTRODUCTION 


The  body-wave  magnitude  of  an  earthquake  was  defined  by 
Gutenberg  and  Richter  as 

mb  =  1oSiq  A/T  +  Q  -3 


where 

A  is  the  P-wave  amplitude  in  millimicrons 

T  is  the  period  in  seconds 

Q  is  the  Gutenberg-Richter  (1956)  distance-depth 
correction. 

This  magnitude  of  course,  depends  on  the  amplitude  re¬ 
corded  at  a  particular  station  observing  the  earthquake.  For 
want  of  a  better  method,  the  arithmetic  mean  of  individual 
station  magnitudes  was  chosen  to  be  representative  of  the 
true  body-wave  magnitude.  This  would  be  correct  if  the  source 
were  assumed  to  be  purely  compressional  (azimuthally  uniform). 

However,  our  knowledge  of  source  mechanisms  has  increased 
to  a  point  where  it  is  now  known  that  the  mechanism  of  most 
earthquakes  may  be  well  represented  by  a  double  couple  source 
mechanism,  Stauder  and  Bollinger  (1964).  This,  in  turn, 
implies  that  the  body-wave  radiation  pattern  is  not  azimuthally 
uniform.  If  we  want  to  be  still  more  accurate,  we  can  also 
include  the  effect  of  a  moving  source.  In  any  case,  the 
arithmetic  mean  of  the  observations  does  not  represent  the 
mean  of  a  radiation  pattern  caused  by  a  fixed  or  moving 
double  couple.  Hence  it  is  desirable  to  find  a  better  way 
of  defining  the  "true"  body-wave  magnitude. 

We  propose  that  a  better  measure  of  the  ’'true"  body- 
wave  magnitude  may  be  obtained  by  taking  the  amplitude  of  a 
purely  compressional  source  whose  radiation  pattern  has  the 
same  area  as  the  observed  radiation  pattern. 


Naturally ,  this  means  that  we  have  to  know  the  source 
parameters,  i.e.,  the  dip  direction,  dip  angle  and  slip  angle 
of  the  fault  plane  and  the  auxiliary  plane. 

In  order  tc  find  the  source  parameters,  two  totally 
.different  methods  have  been  used  to  date.  The  first  is  a 
purely  geometrical  approach,  using  only  the  directions  of 
the  first  motion  of  P-waves.  More  recently  the  study  of  the 
polarization  angle  of  the  S-wave  has  further  helped  to  improve 
estimation  of  the  source  parameters.  However,  as  Stauder  and 
Nuttli  (1S65 )  point  out: 

"Examination  of  multiple  solutions,  i.e.,  by 
several  different  authors  for  one  and  the  same 
earthquake,  as  also  a  comparison  of  groups  of 
fault  plane  solutions  for  a  given  region  evi- 
dense  poor  agreement  and  systematic  differences 
in  solutions  by  different  authors.  This  evi¬ 
dence  requires  that  caution  should  be  exercised 
in  drawing  conclusions  from  the  existing  accumu¬ 
lation  of  published  fault  plane  solutions." 

Although  some  of  the  directions  of  the  first  motions  of 
the  P-waves  are  reported,  no  data  are  readily  available  on 
polarization  angles, and  hence  a  detailed  visual  study  of  the 

seismograms  is  necessary  for  the  analysis  of  source  parameters 
using  the  above  method. 

The  second  method  which  has  been  used  is  spectral  equal¬ 
ization.  This  is  applicable  to  both  body  and  surface  waves 
and  source  parameters  may  be  obtained,  but  with  this  method 
it  is  necessary  to  use  digitized  seismograms  of  a  worldwide 
network,  which  are  not  readily  available. 

The  method  which  we  propose  here  consists  of  using  data 
which  are  readily  available,  i.e.,  reported  first  motion 
directions  together  with  observed  body-wave  magnitudes  as 
giver,  by  the  CSGS  Earthquake  Data  Reports  and  the  ISC  Bulletin. 
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Of  course,  the  method  is  bound  to  be  inaccurate  if  there  are 
large  gaps  in  the  azimuthal  coverage,  but  in  this  case  the 
other  methods  also  fail.  It  is  hoped  that  for  earthquakes 
of  magnitude  5  and  larger,  an  approximate  solution  of  vhe 
source  parameters  may  be  obtained  and  hence  a  better  value 
of  the  mean  body-wave  magnitude  may  be  established. 
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METHOD 


The  method  is  based  on  a  number  of  assumptions: 

1.  Using  only  WWNSS  stations,  it  is  not  necessary  to 
make  instrumental  corrections  since  the  ra^  is  computed  on  the 
basis  of  A/T  and  the  instrument  responses  are  nearly  identical. 

2.  Geometrical  spreading  is  accounted  for  by  the 
Gutenberg-Richter  correction  for  distance  and  depth  of  focus. 

3.  The  major  contribution  to  the  difference  in  the 
observed  variation  of  m^  with  azimuth  is  caused  by  the  source 
radiation  pattern. 

4.  A  double-couple  source  accounts  for  the  radiation 
pattern. 

5.  Sufficient  azimuthal  coverage  exists  in  the  obser¬ 
vations  to  allow  reasonably  unambiguous  determination  of  the 
source  mechanism. 

This  method  is  not  intended  to  replace  the  body-wave 

equalization  procedure  as  given  by  Ben  Menahem,  et  al  (1965). 

It  is  intended  to  provide  a  small  correction  to  the  mean  m,  . 

b 

The  model  taken  is  the  linear  regression  with  zero 
intercept 

B.  =  KA.  +  e. 
l  ii 

M. 

where  B^  is  the  observed  amplitude  (B^  =  10  1) 

K  is  a  constant  to  be  determined 

A^  is  the  calculated  amplitude  (as  given  by  Ben 
Menahem  et  al  (1965) 

e  is  the  error  in  the  ith  observation. 


The  calculated  amplitude  at  the  i’th  station  is  related 
to  the  source  parameters  by: 

Ai  =  a0  +  sin  (0  -  6q)  +  a2  sin  2  (6  -  0q)  + 

b1  cos  (6  -  eQ)  +  b?  cos  2  (6  -  0q) 

where 

a  =  1/4  sin  X  sin  26  (3  cos2iv  -  1) 
o  n 

a^  =  1/2  sin  X  cos  26  sin  2i^ 
b^  =  1/2  cos  X  cos  6  sin  2i^ 
a2  =  -1/2  cos  X  sin  6  sin  2i^ 
b2  =  1/4  sin  X  sin  26  sin  2i^ 

6  is  the  azimuth  (east  of  north)  of  the  station 
as  seen  from  the  epicenter 

0Q  is  the  azimuth  (east  of  north)  of  the  strike 
direction  of  the  assumed  fault  plane. 

6  is  the  dip  angle 
X  is  the  slip  angle 
i^  is  the  take-off  angle 

Figure  22  shows  the  geometry  of  the  source  parameters. 

We  wish  to  use  a  least  squares  procedure  to  determine 
the  source  parameters.  To  do  this,  we  minimize  the  sum  of 

A 

the  squares  of  the  errors,  i.e.,  we  choose  K  so  that 

2 

E  =  Iei  =  “  KA^)  is  minimized. 

Hence  the  estimated  value  for  K  is  given  by 

A 

K  = 
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Since  A^CA)  =  -A^A  +  1800  )  we  determine  which  half-plane 
has  the  smaller  error #  by  choosing 


A 

K 


l?BiAil 

IAi2 


and  then  e .  = 
1 


B.  -  KA. 
i  1 


B.  +  KA. 

i  l 


for  0°  <  A  <  180° 


180°  <  A  <  360° 


Only  stations^for  which  we  have  values  of  m^  were  used  in 
the  estimation  of  K.  In  order  to  include  also  stations  for 
which  mb  is  not  given,  but  the  first  motion  direction  is  given, 
we  assumed  an  observed  value  whose  absolute  value  was  that  of 

the  calculated  value  and  whose  sign  was  the  sign  of  the  observed 
first  motion,  i.e. , 

Bi  =  sign  (STAi)  K | A^ | 

so  that  the  contribution  e^  =  0  if  the  signs  of  the  observed 
and  calculated  values  were  the  same  and  2K|A^|  if  the  signs 
were  opposite.  This  helps  to  insure  that  the  sign  of  the 
larger  amplitudes  will  cause  a  larger  error,  if  they  are 
different  from  the  observed. 

In  order  to  calculate  the  mean  m^,  it  is  first  necessary 
to  define  a  mean  amplitude  for  a  radiation  pattern.  We  intend 
to  use  as  a  mean  amplitude  the  amplitude  of  a  purely  c^mpres- 
sional  (azimuthally  uniform)  source  whose  radiation  pattern 
has  the  same  area  as  the  observed  radiation  pattern: 


-  2 
tta: 
1 


=  1/2  / 


27T  A?  (0)d0. 

O  1 
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I 


t 


This  leads  to 


Aa  =  HrCao  +  1/2  (al  +  a2  +  b!  +  t>2>;1}:L/2 

Now  since  the  coefficients  a  ,  a, ,  a„,  b,  and  are 

0  12  1  2 

functions  of  the  take-off  angle  ih  (which  is  different  for 
each  station),  we  simply  compute  an  arithmetic  mean  of  the 
means  A^,  so  that  finally  the  mean  amplitude  is  given  by 

y  N 

a  =  -4-  l  s. 

~  iii  1 

and  hence  the  mean  m^  is  given  by 

M  =  log1Q 

We  have  thus  calculated  a  new  mean  value  of  m.  which  has  teen 

b 

corrected  for  the  radiation  pattern.  The  standard  error  of 
estimation  of  the  corrected  mb  can  be  calculated  from  the 
t-distribution.  Arbitrarily  taking  our  confidence  limits  at 
90%,  we  find 


K 


.95  w 

U  TT7T 

<1  A2) 

i  =  l  1 
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where  S 


1 

N-l 


l  (B.  -  KA.) 
i  =  1  1  1 


1 

N-l 


2 


N-l 

t  gg  is  the  value  of  the  student's  t-distribution  at 
the  95%  confidence  limit,  for  n-l  degrees  of  freedom. 

Hence  the  upper  and  lower  confidence  limit  for  the  body 
wave  magnitude  is  found  from 


lo'?10  A+ 


K+  N 

where  A.  =  -n-  T  A. 

+  “  i  =  l  1 


K  N 


These  values  are  shown  in  Table  1. 

Algorithm  for  Obtaining  Source  Parameters 

The  initial  strike  angle,  step  size  and  number  of  steps 
are  given  as  input  parameters.  A  coarse  mesh  search  is  carried 
out  to  find  the  best  strike  angle  and  dip  and  slip  angles.  The 
step  size  for  the  dip  and  slip  angles  is  taken  to  be  10°. 

A  fine  mesh  search  is  then  carried  out  keeping  the  strike 
angle  found  in  the  first  statoe  fixed  and  varying  the  dip  and 
slip  angles  to  find  the  least  squares  fit,  in  steps  of  1°. 
Details  are  given  in  the  flow  chart  in  the  appendix  to  this 
report. 
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AMBIGUITY  OF  RESULTS 


When  one  finds  a  fault-plane  solution  using  first  motions 
only,  there  exists  an  ambiguity  in  the  result,  in  that  it  is 
not  possible  to  decide  which  is  the  fault  plane  and  which  is 
the  auxiliary  plane.  Similarly,  using  the  radiation  pattern, 
Ben  Menahem  et  al. ,  (1965)  found  the  amplitudes  and  directions 
of  compressional ,  SV  and  SH  radiation  to  be: 


|(a-S)  (n*S)  5 


<a4>  +  (a*5) 


5 


where 

a  is  the  vector  in  the  direction  of  motion 


n  is  the  vector  normal  to  the  fault  plane. 

Since  the  vector  a,  the  direction  of  motion,  is  normal 
to  the  auxiliary  plane  and  n  is  the  vector  normal  to  the  fault 
plane,  it  can  be  seen  that  the  radiation  pattern  is  unchanged 
by  the  interchange  of  a  and  n;  thus  the  same  ambiguity  which 
exists  for  the  fault-plane  solution,  also  exists  for  the 
radiation  pattern  method. 


Second  Solution  of  the  Fault  Plane 

In  the  previous  section  we  explained  that  in  finding  a 
solution  of  the  source  parameters,  there  exists  another 
solution,  whose  radiation  pattern  is  identical  with  that 
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already  found.  Hence,  since  we  do  not  know  which  is  the  fault 
plane  and  which  is  the  auxiliary  plane,  it  is  necessary  to 
compute  also  the  second  solution.  By  comparing  the  mechanisms 
of  groups  of  earthquakes  from  the  same  area  having  similar  depth 
of  focus,  we  might  be  able  to  see  that  one  of  the  planes  is 
similar  in  all  cases  in  order  to  decide  that  this  is  the  fault 
plane,  assuming  that  the  source  mechanism  does  not  change 
greatly. 


It  is  assumed  that  one  solution 

of  0  ,  6  and  X  has  been  found 
o 

already.  In  Ben  Menahem  et  al 
(1965)  the  three  coordinate  axes 
are  taken  in  the  strike  direc¬ 
tion  (x1>,  a  direction  opposite  to 
the  horizontal  projection  of  the 
dip  direction  (x2)  and  vertically 
upwards  (x^). 

Making  a  transformation  so  that  we 
have  x^  pointing  North,  x'  West  and  x^  vertically  upwards  for  a 
righthand  system,  we  can  describe  the  transformation  of  the 
unit  vectors  by  the  matrix  relation: 


-+• 

“ 

„ 

e. 

cos  0  sin  0  0 

e. 

1 

o  o 

1 

-4*  ” 

e« 

— 

-sin  0  cos  0  0 

-*•  ^ 

2 

o  o 

2 

-*■ 

e3 

—  M 

i 

rH 

O 

O 

_ 1 

e3 

where  0Q  is  the  strike  angle  measured  clockwise  from  North. 

Hence  we  can  describe  the  direction  of  motion,  the  "null" 
vector  and  the  vector  perpendicular  to  the  direction  of  motion 
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»♦ 


by  the  matrix  relation 


a  coa  1  coa#  -aln  »  coa  I  ain»_  coa  X  aln#  »aln  *  coa  4  coat 

o  o  o  o 

(  *  ala  1  coat  -coa  1  coa  I  alnt  ala  1  tint  ♦  coa  1  coa  I  coat 

o  o 

ala  4  aln*.,  -ala  4  coat 

O  W 


The  vector  a  lies  in  the  fault  plane  in  the  direction  of 
motion  and  hence  its  direction  cosines  are  the  direction  cosines 
of  the  auxiliary  plane,  which  is  normal  to  the  fault  plane.  The 
vector  n  is  normal  to  the  fault  plane  and  hence  its  direction 
.cosines  are  the  direction  cosines  of  the  fault  plane.  The  vector 
6  lies  in  the  fault  plane  normal  to  the  direction  of  motion  and 
hence  it  is  the  "null"  vector,  the  intersection  of  the  fault 


ala  1  aln  4  i| 

coa  1  aln  4  .j  (2) 

ooa  4  aj 
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plane  and  the  auxiliary  plane. 

These  relations  are  needed  to  compute  the  slip  angles, 
assuming  that  each  plane  in  turn  represents  the  fault  plane. 

Given  one  solution  of  the  source  parameters  0^,  6^,  A^  we 
wish  to  find  the  second  solution  e2»  52>  *2'  We  can  do  this  by 
remembering  that  a.^  and  n^  of  one  solution  become  n2  and  a2  of 
the  other  solution  and  hence  by  equating  direction  cosines,  we 
have 

cos  A1  cos  ei  +  sin  X1  cos  6±  sin  e±  =  -  sin  &2  sin  02  (3) 

cos  A.^  sin  0.^  -  sin  A.^  cos  6^  cos  61  =  sin  62  cos  (4) 

-  sin  A^  sin  6^  =  -  cos  &2  (5) 

From  these  we  can  derive  the  expression  for  the  constraint 

of  perpendicularity  between  the  two  planes 

tan  61  tan  &2  cos  ( 0 1  -  02>  =  -  1  (6) 

and  the  reciprocal  relations 
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Hence,  having  found  one  solution  for  0^,  6^,  A^,  we 
find  02,  62,  X2  by  finding  6^  from  equation  (8)  where  <$2  is 
defined  to  be  between  0°  and  180°.  In  order  to  find  A  we 
use  equations  (7),  (9),  and  (10)  and  then  we  can  fit  A2  into 
the  correct  quadrant.  Finally,  we  compute  02  by  using 
equations  (6)  and  (9)  and  thus  we  can  fit  02  into  the  correct 
quadrant.  When,  90°  <  62  1  180°,  it  signifies  that  the  dip 
direction  lags  90°  behind  the  strike  direction,  instead  of 
leading  by  90°,  according  to  the  sign  convention  used  here. 
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ESTABLISHING  THE  METHOD 


In  order  to  show  that  the  method  yields  significant 
results,  the  following  test  cases  were  tried: 

1.  Banda  Sea  Earthquake  of  March  21,  1964  (reported  by 
Teng  and  Ben  Menahem  (1965).  For  this  case  the  calculated 
amplitudes  were  also  input  to  check  that  the  method  gives 
back  the  exact  input  values. 

2.  Rat  Island  Earthquake  of  5  February  1965  (origin 
time  09  32  9.3).  The  Rat  Island  earthquake  mechanisms  were 
reported  by  Stauder  (1963).  The  information  for  the  first 
motion  directions  for  the  Rat  Island  sequence  was  kindly 
supplied  by  the  Rev.  William  Stauder,  S.J.  in  a  private 
communication. 

3.  Rat  Island  Earthquake  of  1  October  1965 

4.  Rat  Island  Earthquake  of  15  May  1966. 

5.  Rat  Island  Earthquake  of  4  July  1966. 

6.  Rat  Island  Earthquake  of  22  November  1965. 

7.  Hindu- Kush  Earthquake  of  28  January  1964.  Fault 
plane  solutions  were  reported  by  Hedayati  and  Hirasawa  (1966), 
and  also  by  Ritsema  (1966). 

8.  Alaskan  Earthquake  of  28  March  1964  (origin  time 
03  36  14.2).  The  input  data  was  taken  from  the  Bulletin  of 

the  International  Seismological  Centre.  Solutions  were  reported 
by  Stauder  et  al  (1966)  ,  and  Harding  et  al  (1968). 

9.  Niigata  Earthquake  of  16  June,  1964.  A  fault  plane 
solution  was  reported  by  Hirasawa  (1965). 
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INPUT  DATA 


Input  data  were  collected  mainly  from  the  C  £  GS  Earthquake 
Data  Reports  and  the  Bulletin  of  the  International  Seismological 
Center,  Edinburgh. 

The  first  motion  data  for  the  Banda  Sea  earthquake  were 
taken  from  Teng  and  Ben  Menahem  (1965). 

The  first  motion  data  for  the  Rat  Island  earthquakes  were 
k-ipplifid  by  the  Rev.  William  Stauder,  S.J.  ,  in  a  private 
communication. 

The  first  motion  data  for  the  Hindu-Kush  earthquake  were 
taken  from  Hedayati  and  Hirasawa  (1966). 

The  additional  first  motion  data  for  the  Alaskan  earthquake 
were  supplied  by  Harding  in  a  preprint. 

The  first  motion  data  for  the  Niigata  earthquake  were 
taken  from  Hirasawa  (1965). 

Analysis  of  the  Data 

1.  Banda  Sea  Earthquake  of  the  21st  of  March  1964.  For 
this  earthquake  37  stations  were  used,  12  of  them  with  magnitude 
and  25  with  only  first  motions.  Figure  1  shows  the  observed 
radiation  pattern  and  Table  3  shows  the  input  data.  Figure  2 
shows  the  calculated  radiation  pattern  and  Table  M  gives  the 
results  in  tabular  form.  In  our  final  solution,  36  stations 
had  the  same  calculated  and  observed  first  motion  direction 
and  the  one  disagreement  was  near  a  nodal  value.  The  new 
mean  value  of  m^  came  out  to  be  5.85  as  against  the  CSGS  value 
of  5.8.  The  stars  in  Figure  2  are  calculated  amplitudes  for 
those  stations  which  had  observed  values  and  the  dots  are 
calculated  amplitudes  for  those  stations  which  only  reported 
first  motion  directions. 


# 
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2.  Rat  Island  Earthquake  of  the  5th  of  February  1965  (9h.) 
Figure  3  shows  the  observed  radiation  pattern.  For  this  case 

72  stations  were  used,  13  of  them  with  magnitude  and  59  with 
first  motions  only.  The  results  are  given  in  Figure  4  and 
Table  5.  The  resulting  score  was  72-12,  i.e.,  12  stations 
had  the  opposite  calculated  first  motion.  However,  this  earth¬ 
quake  showed  clearly  that  the  CSGS  value  of  mv  was  too  low 

D 

since  most  of  the  observations  were  near  nodes.  The  new  value 
of  mb  was  6.3  as  against  the  CSGS  value  of  5.9. 

3.  Rat  Island  Earthquake  of  the  1st  of  October  1965.  This 
was  a  case  where  all  the  observed  first  motion  consisted  of 
dilatations .  Figure  5  shows  the  observed  radiation  pattern. 

The  results  are  given  in  Figure  6  and  Table  6.  The  score  was 
80-0.  Although  the  strike  angle  is  quite  far  from  that  found 
by  Stauder  (1968),  the  pattern  of  signs  of  the  first  motion 
are  completely  similar  to  Stauder* s  and  the  solution  was 
determined  primarily  by  the  observed  magnitudes. 

4.  Rat  Island  Earthquake  of  the  15th  of  May  1966.  All 

the  observed  first  motion  consisted  of  compressions.  Figure  7 

shows  the  observed  radiation  pattern.  The  calculated  radiation 

pattern  was  also  governed  mainly  by  the  magnitudes  in  this 

case,  since  the  strike  angle  could  be  varied  considerably 

without  changing  any  of  the  signs  of  the  first  motion.  84 

stations  were  used*  of  which  21  reported  magnitudes  and  the 

rest  only  first  motions.  The  results  are  given  in  Figure  8  and 

Table  7.  This  was  another  case  where  many  of  the  stations 

were  close  to  nodes  and  hence  the  new  value  of  >nL  =  5.96  as 

b 

against  the  CSGS  value  of  5.8  reflects  this. 

5.  Rat  Island  Earthquake  of  the  4th  of  July  1  66.  Figure  9 
shows  the  observed  radiation  pattern.  This  is  a  more  compli¬ 
cated  source  mechanism  and  hence  the  final  score  80-18  shows 
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tf a+  the  fit  was  not  too  good.  22  stations  reported  magnitudes 
and  the  rest  reported  first  motions  only.  One  of  the  diffi¬ 
culties  was  that  no  magnitudes  were  reported  in  the  second 
quadrant,  thus  missing  an  entire  lobe  of  the  pattern.  The 
results  are  shown  in  Figure  10  and  Table  8. 

6.  Rat  Island  Earthquake  of  the  22nd  of  November  1965. 

This  was  another  case  of  compressions  only,  with  no  stations 
reporting  magnitudes  in  the  second  and  third  quadrant.  Figure  11 
shows  the  observed  radiation  pattern.  A  change  of  30°  in  the 
strike  angle  did  not  change  any  of  the  first  motion  directions. 
The  results  are  given  in  Figure  12  and  Table  9. 

Another  run  was  made  using  additional  magnitudes  supplied 
by  D.  Lambert  of  SDL  in  a  private  communication.  The  new 
observed  radiation  pattern  is  given  in  Figure  13.  For  this 
case,  there  were  magnitudes  observed  in  the  third  quadrant, 
but  not  in  the  second  quadrant.  The  strike  direction  did  not 
change,  but  the  dip  and  slip  angles,  and  hence  the  second 
solution,  did  change  somewhat.  The  results  are  given  in 
Figure  14  and  Table  10. 

7.  Hindu-Kush  Earthquake  of  the  28th  of  January  1964. 

This  was  a  complicated  source  mechanism.  The  observed  radiation 
pattern  is  given  in  Figure  15.  Values  of  Logl0  A/T  were 

taken  from  the  ISC  bulletin  and  for  the  first  motion,  only 
those  stations  used  by  Hirasawa  and  Hedayati  (1966)  were 
included.  The  solution  favours  the  one  obtained  by 
Ritsema  (1966).  The  final  score  of  53-14  reflects  the 
inadequacy  of  the  solution.  Again,  one  of  the  difficulties 
is  that  there  are  no  magnitudes  observed  between  115°  and 
295°  azimuth.  The  results  are  shown  in  Figure  16  and 
Table  11. 


-17- 


8.  Alaskan  Earthquake  of  the  28th  of  March  1964  (main 
shock).  The  22  values  of  magnitudf  for  this  event  were 
taken  from  the  ISC  bulletin  and  signs  of  first  motion  were 
taken  from  the  bulletin  and  the  C6GS  Earthquake  Data  Reports. 
The  observed  radiation  pattern  is  given  in  Figure  17.  The 
final  score  obtained,  86-22,  suggests  that  quite  a  few  of 
the  station^  reported  erroneous  first  motions,  since  the 
solution  agress  quite  well  with  that  of  Stauder  and 
Bollinger  (1866).  The  results  are  given  in  Figure  18  and 
Table  12. 

Another  run  was  made  using  additional  first  motions 
supplied  by  Harding  in  a  preprint  of  the  C6GS  on  the  Prince 
William  Sound  Earthquakes  (1968).  The  solution  was  improved 
considerably  and  the  final  score  was  134-17.  The  results 
are  shown  in  Figure  19  and  Table  13. 

9.  The  Niigata  Earthquake  of  the  16th  of  June  1964. 

This  was  another  case  where  m.  given  by  the  C6GS  was  too  low 
and  the  new  value  of  6.5  as  against  6.1  is  correcting  in  the 
right  direction.  The  observed  radiation  pattern  is  given  in 
Figure  20.  The  results,  shown  in  Figure  21  and  Table  14,  do 
not  agree  very  well  with  the  solution  found  by  Hirasawa  (1965). 
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CONCLUSIONS 


Comparison  of  the  fault  plane  solutions  obtained  by  the 
new  method  and  those  obtained  using  first  motions  only  show 
that  indeed  it  is  possible  to  find  reasonable  agreement.  The 
solution  of  the  Hindu- Kush  earthquake  parameters  show  that 
although  better  agreement  with  the  first  motion  directions 
was  obtained  by  Hedayati  and  Hirasawa  (1966),  the  amplitude 
radiation  pattern  does  not  agree  as  well  with  the  observed 
pattern  as  that  obtained  by  using  a  combination  of  first 
motion  directions  together  with  observed  body-wave  magnitudes. 

The  results  of  the  Rat  Island  earthquake  of  the  5th  of 
February,  1965  show  that  the  proposed  method  compensates 
correctly  where  too  many  of  the  observed  magnitudes  were 
near  to  nodes.  This  shows  that  the  proposed  definition  of 
magnitude  is  superior  to  taking  the  arithmetic  mean  of  the 
observations  of  m^. 

The  fault  plane  solution  found  by  this  method  is 
probably  not  as  accurate  as  that  obtained  by  S-wave  data. 
However,  the  new  method  gives  two  checks  on  the  accuracy 
obtained,  the  first  by  simply  the  number  of  differences 
in  sign  between  calculated  and  observed,  i.e.,  the 
"score".  The  second  is  by  examining  the  confidence  limits 
of  the  magnitude. 

The  results  also  indicate  that  in  order  to  take  the 
radiation  pattern  into  account  when  finding  the  "true" 
value  of  m^,  it  is  sufficient  to  use  an  approximate  solution 
to  the  source  parameters.  This  is  possible  by  using  readilv 
available  magnitude  and  first  motion  data  without  laboriously 
having  to  re-examine  records. 

Suggestions  for  Future  Research 

1.  In  order  to  make  best  use  of  the  method  outlined  above, 
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there  still  remains  to  be  determined  xhe  minimum  number  of 
input  data  which  must  be  used  in  order  to  get  a  good  solution. 

For  the  cases  tried  here,  rather  more  points  were  used  than 
really  necessary. 

2.  Instead  of  the  equal  weights  least  squares  solution, 
it  might  be  possible  to  get  some  improvement  by  using  a 
weighted  least  squares  procedure. 

3.  If  more  emphasis  is  to  be  placed  on  the  accuracy  of 
the  fault  plane  solution  obtained  using  the  new  method  discussed 
above,  it  is  possible  to  add  the  S-wave  data  to  the  least-squares 
procedure . 

4.  A  still  better  way  to  define  the  "magnitude"  of  an 
event  may  be  to  use  the  idea  of  moment,  e.g.,  Aki  (1966). 
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Table  3.  Example  of  Input  Data  for  the  Banda  Sea  Earthquake 
of  21  March  1964. 
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Table  4.  Results  for  the  Banda  Sea  Earthquake  of 
21  March  1904. 
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Results  for  the  Rat  Island  Earthquake  of 
5  February  1964. 
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Results  for  the  Rat  Island  Earthquake  of 
1  October  1965. 
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Results  for  the  Rat  Island  Earthquake  of 
15  May  1966. 
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Table  8.  Results  for  the  Rat  Island  Earthquake  of 
4  July  1966. 
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Table  9.  Results  for  the  Rat  Island  Earthquake  of 
22  November  1965. 
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Table  10.  Results  for  the  Rat  Island  Earthquake  of 
22  November  1965  (additional  magnitudes) 
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Table  11. 


Results  for  the  Hindu  Kush  Earthquake  of 
28  January  1964. 
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Table  12. 


Results  for  the  Alaskan  Earthquake  of 
28  March  1964. 
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Table  13. 


Results  for  the  Alaskan  Earthquake  of 
28  March  1964  (additional  first  motions) 
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Results  for  the  Niigata  Earthquake  of 
16  June  1964. 
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Figure  5.  Observed  Radiation  Pattern  for  the  Rat  Island 
Earthquakd  of  1  October  1965. 
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Figure  7,  Observed  Radiation  Pattern  for  the  Rat  Island 
Earthquake  of  15  May  1966. 
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Figure  9.  Observed  Radiation  Pattern  for  the  Rat  Island 
Earthquake  of  4  July  1966. 
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Figure  16.  Calculated  Radiation  Pattern  for  the  Hindu 
Kush  Earthquake  of  28  January  1964. 
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Figure  20.  Observed  Radiation  Pattern  for  the  Niigata 
Earthquake  of  16  June  1964. 
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APPENDIX 


Program  RADPAT  -  Input  Card  Formats  and  Parameters 
Part  I  (only  read  in  once)  Ritsema  i  -  L  table. 

1st  card 

NTAb  Col.  1-5  Number  of  values  in  table  (15) 

2nd  and  subsequent  NTAB/4  cards  (8F10.2) 

ANGCI )  Col.  1-10  t°  (F10.2) 

XIO(I)  11-20  i*  (F10.2) 

Repeat  for  columns  21-40,  41-60,  61-80  as  necessary 
Part  II  (only  read  in  once)  -  Gutenberg-Richter  Distance-Depth 
Correction  Table  (for  use  with  input  values  of  log  A/T). 

X  V 

108  cards,  each  has  format  (4X,  17F4.1),  hence  every  card 
contains  17  fields,  for  a  depth  of  0(25)100,  150(50)700  kilometers. 
Part  III  -  Earthquake  data 

Number  of  earthquakes  to  be  treated  (15) 

Earthquake  identification  (10A8) 

Number  of  stations  (15) 

.LE.O  magnitudes  are  given  as  input  (15) 
.GT.O  log^A/T  is  input  and  the 

Gutenberg-Richter  correction  is 
to  be  applied. 


1st  card 

N  Col.  1-5 

2nd  card 

ID  Col.  1-80 

3rd  card 


M 


ITAB 


Col.  1-5 


6-10 


DEPTH 

11-20 

Depth  of 

the  hypocenter 

(F10.1) 

STRIKEO 

21-30 

Initial  strike  angle  for 

the  search  (F10.1) 

DELSTRI 

31-40 

Increment 

in  sti ike  angle  (F10.1) 

NSTRI 

41-45 

Number  of 

increments  to 

be  used  (15) 

Subsequent  M  cards  contain  the  station  data 


STA(I) 

1-5 

DIST(I) 

11-20 

THETA(I) 

21-30 

MAG (I) 

31-40 

ISIGN(I) 

41-45 

Station  identification  (A5) 

Distance  in  degrees  (F1Q.1) 

Azimuth-epicenter- to-station  in  degrees  (F10. 
Magnitude  or  log1Q  (A/T)  (F10.1) 

First  motion  direction  (15) 

-1  indicates  dilatation 
1  indicates  compression 

0  or  blank  -  if  unknown,  will  be  taken  as 
compression 


COMPUTE  «b 
OS  IBS  GUTENBERG 


END 
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